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Abstract: Understanding solute transport behaviors of deep soil profile in the Loess Plateau is helpful for 
ecological construction and agricultural production improvement. In this study, solute transport processes 
of a deep soil profile were measured by a conservative tracer experiment using 25 undisturbed soil cores 
(20 cm long and 7 cm diameter for each) continuously sampled from the surface downward to the depth 
of 500 cm in the Loess Plateau of China. The solute transport breakthrough curves (BTCs) were analyzed 
in terms of the convection-dispersion equation (CDE) and the mobile-immobile model (MIM). Average 
pore-water velocity and dispersion coefficient (or effective dispersion coefficient) were calculated using 
the CDE and MIM. Basic soil properties and water infiltration parameters were also determined to 
explore their influence on the solute transport parameters. Both pore-water velocity and dispersion 
coefficient (or effective dispersion coefficient) generally decreased with increasing depth, and the 
dispersivity fluctuated along the soil profile. There was a good linear correlation between log-transformed 
pore-water velocity and dispersion coefficient, with a slope of about 1.0 and an average dispersivity of 
0.25 for the entire soil profile. Generally speaking, the soil was more homogeneous along the soil profile. 
Our results also show that hydrodynamic dispersion is the dominant mechanism of solute transport of 
loess soils in the study area. 


Keywords: solute transport; loess soil; pore-water velocity; dispersion coefficient; hydraulic conductivity; Loess 
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1 Introduction 


Solute transport is an important mode of nutrient flow to plant roots and agro-chemicals to 
groundwater (Shao et al., 1998; van der Zee and Leijnse, 2013). Understanding the behavior of 
soil water and soluble substances is critical for increasing nutrient use efficiency and protecting 
shallow groundwater systems. However, field-scale soil solute transport is complicated and 
affected by hydraulic properties and heterogeneity of the soil (Huang et al., 1995; Seuntjens et al., 
2002; Mallants, 2014). Thus, both theoretical and experimental models have been developed for 
investigating the behavior of solute transport in soils (e.g., Tsuboyama et al., 1994; Seuntjens et 
al., 2002; Gao et al., 2010; Koestel et al., 2013; Venkatraman et al., 2014). Theoretically-driven 
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mathematical models such as the convection-dispersion equation (CDE) can be used to simulate 
nonreactive solute transport through saturated porous media under laboratory and field conditions 
(Lapidus and Amundson, 1952; Nielsen and Biggar, 1962). Also the mobile-immobile model 
(MIM) can be utilized to describe anomalous transport in heterogeneous porous media (van 
Genuchten and Wierenga, 1976). 

Theoretical models are widely used to conduct experiments on various elements of solute 
transport in soils (Lennartz, 1999; Zhou et al., 2011; Chotpantarat et al., 2012). However, the 
well-controlled boundary conditions required for theoretical models are difficult to establish at the 
field scale. Undisturbed soil columns with well-preserved natural layers can provide useful 
information on the behavior of solute transport at local scales (Seyfried and Rao, 1987; Mallants, 
2014). The depth of most undisturbed soil columns is within few to tenths of centimeters, but the 
depth of some large columns can reach up to =>100 cm (Vanclooster et al., 1995; Garré et al., 2010). 

The Loess Plateau of China has the largest deposit of loess soil (both in area and depth) and 
also has the most severe soil and water loss, resulting in ecosystem degradation (Shi and Shao, 
2000). Thus, efforts have been made to investigate the spatial distribution of soil water and 
nutrient at different scales to restore and reconstruct the degraded ecosystem (Fu et al., 2000; Gao 
et al., 2015; Li et al., 2015). Researches on the processes of solute transport in soils of the Loess 
Plateau has focused mostly on repacked (Zhou et al., 2011) and undisturbed (Jiang et al., 2012; 
Yang et al., 2013) soil columns sampled from the surface downward to the depth of less than 100 
cm. Little has been done on solute transport behavior in deep soil profiles. Some studies focused 
on hydrological properties of deep soil layers indicated that soil hydrological properties varied 
little at the depth of >500 cm (Wang et al., 2013; Li et al., 2016). As solute transport depends 
largely on water movement, evaluating the solute transport processes of soil profile at the depth of 
0-500 cm could have significant implications for the management of soil nutrient and 
groundwater quality. 

In this study, we conducted a solute transport experiment in undisturbed soil cores collected in 
vertical successions from a deep soil profile of the Loess Plateau, and characterized the effects of 
soil properties and water infiltration parameters on solute transport processes. The parameters for 
the infiltration of water and solute transport breakthrough curves (BTCs) were determined under 
steady-state flow conditions. The BTCs were fitted using both the CDE and the MIM. Some basic 
soil properties were used to explain the variations in solute transport processes of the soil profile. 


2 Materials and methods 


2.1 Undisturbed soil core sampling 


The undisturbed soil cores were collected in a corn field at Yangling (34°30'N, 108°07'E) situated 
in the Loess Plateau, China. The annual mean temperature at the study site is 12.8°C and the 
mean annual precipitation is 530 mm. The soil in the field is Lou soil (Eum-orthic Anthrosol 
classified using Chinese Soil Taxonomy; Wang et al., 2006). The undisturbed soil cores were 
collected from 28 September to 2 October in 2017, when the corn was ripe and ready for harvest. 
A soil pit (150 cm long, 100 cm wide and 500 cm deep) was dug and one of the vertical faces was 
smoothened with a shovel. The soil core sampling was done using a transparent acrylic cylinder 
(7 cm in internal diameter and 24 cm in height). The thickness of the cylinder wall was 5 mm with 
a sharp bottom end for cutting through soil profile. The cylinder was gently driven into the soil at 
3-cm depth each time using a rubber hammer and then removed. The soil outside the cylinder was 
scrapped off. This process was repeated until the soil core was about 20 cm into the soil, and the 
depth was determined by sticking a tape on the wall of the cylinder. The soil core was cut off and 
covered in a piece of filter paper, a perforated plate and a lid with a drainage hole at the bottom 
end. The top end was plugged with soil in a plastic bag to prevent the sample from falling off. A 
total of 25 undisturbed soil cores (20 cm long and 7 cm diameter) were collected in this way from 
the surface down to the depth of the smoothened face (i.e., 0-500 cm). Then, the soil cores were 
transported to the laboratory. 
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2.2 Water infiltration and solute transport experiment 


In the laboratory, the soil cores were placed on custom-made shelves. The surface of each soil 
core was covered with a piece of filter paper to prevent the formation of crust during watering. 
Small sampling bottles were put under the soil cores to collect drainage water. Deionized water 
was applied to the top of the cores using Mariotte bottles to maintain a constant water head, which 
was recorded for each case for the latter calculation. The water table in the Mariotte bottle was 
measured to calculate the infiltration rate with short-time intervals for the early period and 
long-time intervals for the later period. The infiltration experiment was terminated when the 
wetting front reached the bottom of the core. Then the ponding water at soil surface was removed 
and 0.2 mol/L NaCl solution was used as tracer to seamlessly replace the deionized water. 
Subsequently, the effluent was collected continuously in a 25 mL graduated cylinder at fixed 
intervals and titrated with 0.1 mol/L AgNO; solution to determine Cl concentration until it was 
close to 0.2 mol/L. 


2.3 Soil property determination 


In order to measure the basic physical properties, we collected another 50 undisturbed soil 
samples from the same profile at 10 cm intervals close to the soil core using a cutting ring with 5 
cm diameter. The saturated hydraulic conductivity (Ks) was determined using the constant-head 
method. Soil samples were oven-dried at 105°C for 48 h after Ks measurement to calculate soil 
water content and bulk density. The disturbed soil samples were also collected at the same depth 
interval (i.e., 10 cm interval) for the determination of particle size distribution by laser diffraction 
method. 


2.4 Model description 


The CDE and MIM were used to simulate the solute transport processes. The one-dimensional 
equilibrium CDE is given as Equation 1 (Lapidus and Amundson, 1952): 
2 
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where R is the retardation factor; C is the resident concentration (g/cm*); t is the time period (min); 
D is the dispersion coefficient (cm?/min); v is the pore-water velocity (cm/min); and x is the 
distance (cm). For non-reactive tracers, R is equal to 1.0. The MIM for non-equilibrium transport 
given by van Genuchten and Wagenet (1989) is expressed as follows (Eqs. 2 and 3): 
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where Qis the volumetric soil water content (%); a is the first-order mass transfer coefficient 
(min™!); and the subscripts m and im refer to the mobile and immobile regions, respectively. The 
measured BTCs were fitted for the CDE and MIM to estimate solute transport parameters using 
the CXTFIT non-linear fitting procedure (Toride et al., 1995). The effective dispersion coefficient 
(Dett) was derived from Equation 4, and it was introduced in the MIM to compare the solute 
transport parameters between CDE and MIM (van Genuchten and Dalton, 1986). 
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where the first term of the right part of Equation 4 denotes the hydrodynamic dispersion in the 
mobile phase, and the second term reflects the source-sink effects of solute transfer into and out 
of the immobile regions. 


2.5 Data analysis 


The solute transport parameters were estimated using the CXTFIT computer program (Toride et 
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al., 1995). The figures were plotted with Sigmaplot 12.5 (Systat Software Inc., UK). 


3 Results and discussion 


3.1 Soil profile properties 


The soil bulk density along the soil profile was in the range of 0.91—1.35 g/cm*, with an average 
of 1.10 g/cm? (Fig. la). The average soil bulk density in the 0-20 cm soil layer was 1.29 g/cm’, 
similar to the value of undisturbed soil in the Loess Plateau reported by Jiao et al. (2011). Soil 
compaction initially decreased and then increased from the top to the bottom of the soil profile. 
Wang et al. (2008) observed bulk density distribution in the 0-100 cm loess soil profile and found 
that the 10-20 cm soil layer (the so-called plow pan) had the highest compaction. 

The vertical distribution of initial gravimetric water content in the soil profile is plotted in 
Figure 1b. The soil moisture decreased with increasing depth for the surface 0-20 cm soil layer 
and then increased with increasing depth for the 20-100 cm soil layer. The soil moisture 
decreased slowly with increasing soil depth below 100 cm depth, with no significant differences 
between individual soil layers except for the 390-400 cm soil layer (which could have been due 
to error). Generally speaking, the shallow soil layer was low in moisture because of the 
evaporation (Heathman et al., 2003; Wang et al., 2008). The high soil moisture in the surface soil 
layer in this study may be due to the reason that soil sampling was done after low precipitation 
events. 

The clay content in soils increased at first and then decreased later (Fig. 1c). The maximum 
clay content was observed at the 30—40 cm soil layer. Soils with high clay content usually have a 
high water retention capacity, suggesting that water content in the surface soil layer was largely 
influenced by climatic and vegetation conditions in the study area. 
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Fig. 1 Vertical distributions of basic soil physical properties along the soil profile 


3.2 Water flow characteristics 


The final infiltration rate of each soil core is shown in Figure 2. No effluent was collected for soil 
core samples from the 360-380 cm soil layer after 30 d of application of water. Therefore, the 
infiltration rate of this soil layer was assumed to be 0. It is shown in Figure 2 that the soil column 
above the 200 cm soil depth had a stronger infiltration capacity than the soil column below the 
200 cm soil depth. Compared with the infiltration rate calculated by Wang et al. (2008), there was 
a markedly high infiltration rate for the 100—200 cm soil layer in this study. The macro-pores in 
this soil layer could explain the high infiltration with local peak values of 60—80 cm soil layer. 
However, Kang et al. (2002) observed a much greater infiltration rate of 0.16 cm/min for the 0—60 
cm soil layer in a winter wheat field in the southern region of the Loess Plateau, and Zhou et al. 
(2010) noted an infiltration rate of 0.12 cm/min for the 0-10 cm soil layer in the northern region 
of the Loess Plateau. The discrepancy may be due to the differences in land use types and 
vegetation types. 
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The vertical variations of K, along the soil profile shown in Figure 3 indicated that K, was low 
both in the 0-200 and 350-500 cm soil layers. The relatively high K, values in the 200-350 cm 
soil layer and at the 480 cm soil depth may be due to the large soil pores. The differences between 
the distributions of final infiltration rate and K, reflected a high spatial variability of soil water 
infiltration in the study area. 

Final infiltration rate (cm/min) 
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Fig.2 Vertical distribution of final infiltration rate along the soil profile 
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Fig.3 Vertical distribution of saturated hydraulic conductivity (K,) along the soil profile 


3.3 Solute transport processes 


The coefficients of determination (R°) for curve fitting the BTCs using both CDE and MIM were 
greater than 0.9, indicating that the fitting results were reliable. Average pore-water velocity is 
independently estimated and plotted in Figure 4a, showing a decreasing trend with increasing soil 
depth in general. The average pore-water velocity was higher in the sampled soil cores at the 
depth of 0-200 cm than at the depth of 200-500 cm. The result was in a good agreement with the 
distribution of final infiltration rate. The pore-water velocity varied dramatically at the 0-200 cm 
soil depth, but it was distributed uniformly at the 200-500 soil depth, except for the last 20 cm 


JOURNAL OF ARID LAND 


soil depth. As reported by Mallants et al. (1994), soils near the surface are more heterogeneous. 

The mobile-phase velocity (vm) was obtained using the MIM (v,=v/(6,/@), where vmis the 
mobile-phase velocity (cm/min), v is the average pore-water velocity (cm/min), @ is the 
volumetric soil water content (%) and the subscript m refers to the mobile region). The 
distribution of mobile-phase velocity along the soil profile is plotted in Figure 4b. The high values 
at the 100, 200, 420 and 460 cm soil depths could be due to the very low mobile water fraction, 
where nearly all the water was immobile. Without taking this into consideration, the mobile-phase 
velocity was less variable than the pore-water velocity. Generally, the mobile-phase velocity 
decreased with increasing soil depth and became almost constant starting from the 220 cm soil 
depth downward to the 360 cm soil depth. 
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Fig. 4 Vertical distributions of average pore-water velocity (v, a) and mobile-phase velocity (vm, b) along the 
soil profile 


Both of the dispersion coefficient and effective dispersion coefficient were generally larger at 
the shallow depth than at the deep depth (Fig. 5). The dispersion coefficient varied considerably at 
the depth of 0-200 cm. However, it became almost constant at the depth of 200-500 cm, except 
for the deepest soil depths. Soils in the upper region of the soil profile generally had higher 
porosity values and looser structures, which facilitated solute transport. Soils were more compact 
at the deep depth, causing high tortuosity of the processes of solute transport. Similarly, the large 
dispersion coefficient values for soil cores in the 480-500 cm soil column were due to the 
presence of macro-pores, which could result in a high water mobility. Generally speaking, the 
dispersion coefficient is proportional to the pore-water velocity, thus an increase in pore-water 
velocity always results in an increase in dispersion coefficient. This consistency was also 
observed in this study, where both pore-water velocity and dispersion coefficient were high at the 
depth of 0-200 cm. 

Comparison of dispersion coefficient pertaining to CDE and MIM (dispersion coefficient and 
effective dispersion coefficient, respectively) showed that the values derived from the former 
were in some cases larger than those derived from the latter. This was expected because solute 
exchange process in MIM was the equivalent of dispersion in CDE (Wang and Shao, 2007; Zhou 
et al., 2011). Therefore, dispersion coefficient fitted by the MIM accounted for only 
hydrodynamic dispersion. However, the effective dispersion coefficient values at the 160 and 240 
cm soil depths were larger than the corresponding dispersion coefficient values. This was due to 
the exchange of solutes between mobile and immobile regions, which could result in relatively far 
solute travel distances, thus requiring a larger effective dispersion coefficient value to describe it. 

Solute transport in structured soils can approach equilibrium CDE when local equilibrium 
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Fig. 5 Vertical distributions of dispersion coefficient (D) and effective dispersion coefficient (Det) along the soil 
profile 


assumptions are met (Parker and Valocchi, 1986). With low pore-water velocity, the mean 
residence time is larger than the characteristic time for solute exchange between mobile and 
immobile phases. Therefore, the validity of local equilibrium assumptions is supposed to improve 
for compact soils. In this study, comparison between dispersion coefficient and effective 
dispersion coefficient indicated a close agreement for soil cores from 240 to 460 cm depths. Soil 
bulk density increased downwards starting from 200 cm depth (Fig. 1), indicating that deep soils 
are more compact than shallow soils. For the soil cores at the deep depth (200-500 cm), solute 
transport could be described by means of CDE model due to the compaction of soil. 

In this study, the estimated values of dispersion coefficient increased with increasing 
pore-water velocity (R?=0.84, P<0.01). For each soil core, dispersivity (A, cm) was used to 
describe the relationship between pore-water velocity and dispersion coefficient (A=D/v, where D 
is the dispersion coefficient (cm?/min) and v is the pore-water velocity (cm/min)). The calculated 
value of dispersivity fluctuated along the soil profile, with a range of 0.16-5.31 cm and an 
average of 1.61 cm. Dispersivity is generally regarded to be characteristic of the entire medium 
(Bear, 1972). It has typical values of 0.1-2.0 cm for repacked homogeneous soil columns, which 
are one or several times higher than those for heterogeneous soil columns (Huang et al., 1995). 

Values of pore-water velocity and dispersion coefficient for all soil cores were combined to 
evaluate the relationship between the two parameters for the whole soil profile. A linear 
regression equation (Eq. 5) was used to estimate the dispersivity (A) and exponent (n). 

D=4v", (5) 
where n is an empirical parameter, usually in the range of 1-2 (Saffman, 1959; Beven et al., 1993). 
Log-transformed pore-water velocity and log-transformed dispersion coefficient plotted in Figure 
6 showed a good linear correlation. The values of dispersivity and exponent were determined 
from the regression correlation between log-transformed pore-water velocity and log-transformed 
dispersion coefficient. By the combination of all the soil cores, ensemble mean values of 
dispersivity and exponent were obtained, with the values of 2.25 cm and 1.1561, respectively. 
Such a small dispersivity value for undisturbed soil column was also compiled by Jury (1985) and 
observed by Mallants (2014). The small dispersivity value probably reflected the relative 
homogeneous nature of loess soil profile, especially at deeper soil layers. The value of exponent 
was close to 1.0, indicating that hydrodynamic dispersion is the main solute transport process in 
the study area. A similar exponent value (n=1.2) was observed by Jaynes et al. (1988) in a 
bromide leachate experiment under continuous flood conditions. Due to the heterogeneity of the 
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soil, the correlation between vertical distributions of solute transport parameters and soil 
properties was not firmly conclusive. However, the above discussion suggested that the spatial 
fluctuations in solute transport parameters can be explained qualitatively by specific soil physical 
properties. 
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In(D)=1.1561In(v)+0.8123 
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8 7 -6 5 + 3 -2 -l 
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Fig. 6 Correlation between log-transformed pore-water velocity (In(v)) and log-transformed dispersion 
coefficient (In(D)) 


4 Conclusions 


A conservative tracer experiment was conducted to characterize the solute transport processes of a 
deep loess soil profile (0-500 cm depth) in the Loess Plateau. Average pore-water velocity was 
determined independently from the experiment and dispersion coefficients were calculated by 
fitting the CDE and MIM to the solute concentration vs. time data. The final infiltration rate 
initially increased and then decreased with increasing soil depth, showing a negative correlation 
with soil bulk density. The BTCs data for all soil cores can be explained by both CDE and MIM. 
The measured pore-water velocity had larger values in the shallow soil layers, with fluctuations 
driven by macro-pores in the soil. Both dispersion coefficient and effective dispersion coefficient 
generally decreased with increasing soil depth, except for few soil layers. The close agreement 
between the two parameters (dispersion coefficient and effective dispersion coefficient) below the 
200 cm soil depth indicated that the deep soil layers were largely homogeneous. The dispersivity 
varied considerably over the entire soil profile and had a magnitude smaller than that for some 
typical soil cores. The log-transformed pore-water velocity and log-transformed dispersion 
coefficient exhibited a linear correlation, and the result indicated that the soil was relatively 
homogeneous. Our results further showed that hydrodynamic dispersion is the main spreading 
mechanism of solute transport of loess soils in the study area. 
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